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The question of how many entangled or, respectively, sep- 
arable states are there in the set of all quantum states is con- 
sidered. We propose a natural measure in the space of density 
matrices q describing iV-dimensional quantum systems. We 
prove that under this measure, the set of separable states pos- 
sesses a nonzero volume. Analytical lower and upper bounds 
of this volume are also derived for N — 2 x 2, and iV = 2 x 3 
cases. Finally, numerical Monte Carlo calculations allow to 
estimate the volume of separable states providing numerical 
evidence that it decreases exponentially with the dimension of 
the composite system. We have also analyzed a conditional 
measure of separability under the condition of fixed purity. 
Our results display a clear dualism between purity and sep- 
arability: entanglement is typical for pure states, while sep- 
arability is connected with quantum mixtures. In particular, 
states of sufficiently low purity are necessarily separable. 

03.65.Bz, 42.50.Dv, 89.70,+c 



I. INTRODUCTION 

The question of quantum inseparability and entangle- 
ment of mixed states has attracted much attention re- 
cently. This problem is, by far, more complicated than 
the analogous one for pure states [El an involves subtle 
effects like "hidden nonlocality" M or "distillation of en- 
tanglement" H[l]]. Generally speaking, one is interested 
in inseparable states as the states containing Einstein- 
Podolsky- Rosen (EPR) correlations. In fact, all insep- 
arable mixed states have nonzero "entanglement of for- 
mation" j^] which means that to build them a nonzero 
amount of pure entangled states is needed. In particular 
if a source emits pairs of particles in unknown pure states, 
so that they form a quantum ensemble described by an 
inseparable density matrix, then it follows that the source 
must emit some entangled pairs with a nonzero probabil- 
ity. In this sense the inseparable mixed states can be 
vieved as entangled in correspondence to the entangled 
pure states. 

One of the fundamental questions concerning those 
subjects is to estimate how many entangled (disentan- 



gled) states exist among all quantum states. More 
precisely, one can consider the problem of quantum 
separability-inseparability from the measure theoretical 
point of view, and ask about relative volumes of both 
sets. There are three main reasons of importance of this 
problem. The first reason, of some philosophical impli- 
cations, may be contained in the questions "Is the world 
more classical or more quantum! Does it contain more 
quantum correlated (entangled) states than classically 
correlated ones?" . The second reason has a more practi- 
cal origin. Analyzing some features of entanglement one 
often has to rely on numerical simulations. It is then 
important to know, to what extend entangled quantum 
states may be considered as typical. Finally, the third 
reason has a physical origin. Physical meaning of sepa- 
rability has been recently associated with the possibility 
of partial time reversal Q (see also 0). Separable states 
of composite systems allow time reversal in one of sub- 
systems, without loosing their physical relevance. How- 
ever, for a system of a dimension N > 8 the fact that a 
state admits partial time reversal is not sufficient to as- 
sure separability, and counterexamples have been found 
||. Moreover, it has been recently shown that none of 
those counterexamples can be distilled to a singlet form 
Therefore, it seems pertinent to investigate how fre- 
quently such peculiar states appear. At first glance it 
seems quite likely that such states form a set of mea- 
sure zero, and that from a measure theoretical point of 
view the set of separable states and the states that admit 
partial time reversal have equal volumes. 

In this paper we make an attempt to answer at least 
the two first of the above formulated questions. We also 
give a qualitative argument why the last conjecture fails. 
To this aim we propose a simple and natural measure on 
the set S of density matrices acting on a finite dimen- 
sional Hilbert space TL. Using this measure we estimate 
the relative volume of the set of separable states S sep . 
The upper (lower) bound on this volume is obviously the 
lower (upper) bound on the relative volume of the set of 
inseparable (entangled) states Si nsep = S\ S aep . 

The paper is organized as follows. Section II contains 
our definition of the natural measure in S. In Section 
III we recall basic definitions of separable states, and 
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prove that for any compound system S, the volume of 
S sep is nonzero regardless the number of subsystems it 
contains and its (finite) dimension. This is achieved by 
proving the existence of a topological lower bound of this 
volume. Better lower bounds are also calculated analyti- 
cally by analyzing the relation between the purity of the 
state and its separability. In Section IV analytic upper 
bounds on the volume of S sep are found. The study of 
inseparable states with positive partial transposition is 
presented in Section V. In Section VI we present the es- 
timates on the volume of separable states obtained by the 
Monte Carlo numerical simulations. This section is self- 
contained, and includes also a simplified corollary of the 
results of Sections II-V, and a discussion of the dualism 
between purity and separability. We conjecture that the 
volume of separable states decreases exponentially with 
the dimension of the Hilbert space. Finally, Section VII 
contains our conclusions and open questions. 

The reader should note the that the Sections II-V have 
a rather formal mathematical character. The results of 
these Sections provide a rigorous base for the numerical 
calculations of Section VI, but the detailed knowledge of 
the proofs is by no means necessary to understand the 
main message of the paper. The reader who is not inter- 
ested in such rigorous proofs of the presented results, may 
well skip Sections III-V, and go straight to the Sections 
VI and VII. 



II. NATURAL MEASURE OF THE SET OF 
QUANTUM STATES 

Let us consider a set of states in a V-dimensional 
Hilbert space TL. In particular, TL may describe a compos- 

m 

ite system with m component subsystems: TL = <£> TLi, 

i=l 

where n£Li N { = N. 

An operator g acting on TL describes a state if Trg = 1 
and if g is a positive operator, i.e. 

Tr(gP) > (1) 

for any projector P. Any state represented by a den- 
sity matrix g can in turn be represented by its spectral 
decomposition: 

N N 

g=J2 A nPn, ^A„ = l, A„>0, (2) 

n—1 n—1 

where P n form a complete set of orthogonal projectors. 
Thus the set of states can be viewed as a Cartesian prod- 
uct of sets: 

S = VxA. (3) 

The set V denotes the family of complete sets of or- 

thonormal projectors {Pi}f-\, Yl n =i ^" = ^> wnere I 
is the identity matrix. There exist the unique, uniform 



measure v on V induced by the Haar measure on the 
group of unitary matrices U(N). Integration over the set 
V amounts thus to the integration of the corresponding 
angles and phases in V-dimensional complex space that 
determine the families of orthonormal projectors (or, al- 
ternatively speaking the unitary matrix that diagonalizes 

The symbol A in Eq. (Q) represents there the set of all 
A„'s, which is a subset of the N—1 dimensional linear 
submanifold of real space 1Z N , defined by the trace condi- 
tion, y^ =1 A„ = 1. Geometrically, A is defined as a con- 
vex hull (i.e. a set of all convex combinations of the edge 
points) A = convjxi e : Xi = (0, L;, 0), i = 
1, N}. Since the simplex A is a subset of the N—1 
dimensional hyperplane, there exist a natural measure on 
A which is defined as a usual normalized Lebesgue mea- 
sure jCjv-i on Ti- 1 ^" 1 ■ More specifically, any measurable 
function /(.) of Ai, . . . , Ajv can be integrated with this 
measure: 

,1 ,1 N 

— dAi... dA A r/(A 1 ,...,A iV )(J(VA n -l) = 

V s Jo JO j 

1 rl rl N-l 

— dA^.. dAjv-i/(Ai,...,Aiv-i,l- V A n ), (4) 

V s Jo JO q 

where the normalization constant V s equals to the vol- 
ume of the set A in TZ 1 ^ -1 , whereas S(.) denotes the 
Dirac's delta distribution. The two above discussed mea- 
sures induce a natural measure on S: 

H = vx C N -\- (5) 



III. VOLUME OF THE SET OF SEPARABLE 
STATES 

A. Preliminaries - separable states 

Throughout this paper we shall assume that the 
Hilbert space of the considered quantum system has an 
arbitrary but finite dimension. To make further consid- 
erations more clear, we start from the following notation 
and definitions. Recall first that the space A of operators 
acting on TL constitute a new Hilbert space (a so called 
Hilbert-Schmidt space) with the scalar product: 

(A,B) = Tr(B^A). (6) 

It induces a natural norm (a trace norm): 

= ^Tr(AU). (7) 

which, according to the condition dimTL < oo is topolog- 
ically equivalent to all other norms on A, in particular to 
the norm \ \A\\' = Tr|A|. Furthermore, let us recall that: 
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Definition 1 The state g acting on the Hilbert space 
TL = Tli <8> Ti-2 is called separable [] if it can be approxi- 
mated in the trace norm by the states of the form 

k 

Q = ^2piQi®Qi (8) 
i=l 

where gt and Qi are states on Ti\ and TL2 respectively. 

Usually one deals with a finite dimensional Hilbert space 
dim7i = N. For this case it has been shown || that 
any separable state can be written as a convex combina- 
tion of finite product pure states, i.e. in those cases the 
"approximation" part of the definition is redundant. 

It has also been shown that the necessary condition 
for separability of the state g is positivity of its partial 
transposition g T2 . The latter is defined in an arbitrary 
orthonormal product basis \fi) (g> \fj) as a matrix with 
elements: 

Qmm',nn> = (/"'I ® (fm> |£> T2 |/n) ® |/„') = Qmn>,nm>- (9) 

Although the matrix g 1 " 2 depends on the used basis, its 
eigenvalues do not. Consequently, for any state the above 
condition can be checked using an arbitrary product or- 
thonormal basisQ. For systems of dimensions 2x2 and 
2x3 the partial transposition condition is also a sufficient 
one and thus the set of separable states is completely 
characterized by this condition. 

The definition of separable states can be easily general- 
ized to systems composed of more than two subsystems: 

Definition 2 The state g acting on the Hilbert space 

m 

Ti = ® TLi is called separable if it can be approximated 
2=1 

in the trace norm by the states of the form 

k 

Em . 
Pi ® Qi (10) 

where g\ are states on Tti . 

Straightforward generalization of the proof about decom- 
position from Ref. || gives us the possibility of omitting 
the approximation part in the definition: 

Lemma 1 Any separable state g of a system composed 
by m subsystems can be written as: 



1 The presented definition of separable states is due to 
Werner jD], who called them classically correlated states. 

2 As the full transposition of a positive operator is also pos- 
itive, positivity of the partial transposition g T2 is equivalent 
to positivity of the partial transposition g Tl (defined in anal- 
ogous way). 



k 

= Y,P iP prod, k<N 2 (11) 
where Pp rod are pure product states having the m- 

m 

decomposable form ® Pi, where Pi are projectors acting 
1=1 

on Tli. 

It is worth mentioning that minimal decompositions with 
k = N can be always found for N — 4 pjl^j . 

B. Existence of nonzero lower bound for the volume 
of separable states 

We shall prove now that the volume of the set of sep- 
arable states is nonzero independently of the dimension 
of the Hilbert space and the number of subsystems m 
composing it. 

For our purposes we first prove the simple lemma: 
Lemma 2 If the hermitian operator A € A satisfies 

TCI TCI 

(A, <X> Pi) = for any product projectors <X> Pi, then 

i—1 i—1 

it is a trivial zero one. 

Proof - Let us consider an arbitrary orthogonal (in the 
sense of scalar product (Q)) product hermitian basis in 
the space of operators A, i.e. a basis such that any of its 
elements is a product of hermitian matrices (for instance, 
in the 2x2 case the basis could consist of products of 
Pauli matrices a n ® er m , with n, m = 0, 1, 2, 3; <7o = / ). 
Hermicity assures that any element of the basis can be 
written as a real combination of product projectors. Any 
coefficient of the expansion of A in this basis is given by 
the scalar product @ of A and the corresponding ba- 
sis element. From the general assumption of vanishing 

m 

of formulas of type (A, ® Pi) = 0, we obtain immedi- 

i—l 

ately that all expansion coefficients must vanish. Hence 
A must be equal to the zero operator. 

We can now propose the following theorem : 

Theorem 1 Let A c be a simplex defined as A e = 
conv{y ; e TZ N : yi = ex ; + (1 - e)z I :i = 1,...,N; zi = 

(j?, jf )}- Let define a set Q e = V x A e . Then there 
exists some positive e such that Q e C S sep , where S sep 
represents the set of separable states. 

The meaning of the above theorem is straightforward. 
It says that all states in the sufficiently small neighbor- 
hood of the maximally mixed state gi = -kl (which 
is represented in A as a point zi = ( , . .. , jr ) for any 
chosen spectral decomposition of unity) are necessar- 
ily separable. Note that by definition, the simplex A e 
has edges e times smaller than A, so that its volume 
/i(A e ) = e 7V_1 /i(A) = e w_1 , since according to our nor- 
malization /i(A) = 1. 

Proof .- Suppose, conversely that for any positive e the 
set Q e contains some inseparable state gi nse p ■ It is easy to 
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see that then there must exist a sequence of inseparable 
states g" nsep convergent to the maximally mixed state 
qi . According to Lemma 1 and Theorem 1 from Ref . ]Tc| ] 
there exist a sequence of operators A n separating the 
states g™ nsep from the state gi in the sense that for any 
n it holds: (A n , Q? nsep ) < and (A n , Qz) > 0. Moreover, 
from the quoted results it follows that (A n , a) > for 
any a € S sep . Let us normalize the operators A n by 
introducing A n = A n /\\A n \\ ( \\A n \\ = y/(A^A)). These 
operators satisfy: 

(A n , g™ nsep ) < 0, (A n , a) > for any a G S sep . (12) 

In particular it holds {A n , qi) > 0. ^From construction 
the sequence A n belongs to the sphere in the finite di- 
mensional space A. As the latter is a compact set, the 
sequence includes some subsequence A n ^) which is con- 
vergent to some nonzero operator A* = 1). ^From 
( |l2| ) and continuity of scalar product it follows that the 
limit operator also satisfies: 

(A*, a) > for any a G S sep . (13) 

Now using (|l2|) and the Schwarz inequality we obtain: 

< (A»(fc), Qi) = (An(k), Qi ~ Qlnslp) + (An(k), otnlp) 



< (A n{k) , qi - g"^l p ) < \ \A n{ k)\\\\gi - g 



n(k) , 
insep I 



\Qi - Q, 



>(*0 



(14) 



Taking the limit with respect to k we obtain: 
Tri 



(A*,qi) = lim {A n(k) , Ql ) = 0. (15) 

k — >oo 



Hence A* is traceless, which is in contradiction with 
(|l3|). Indeed, if the operator A* is to be nontriv- 
ial (the construction implies its unit norm) then there 

m 

must exist some product state P pro d = <8> P % such that 

i=l 

(A*, Pprod) (see lemma 2). Since, on the other hand, 
one requires the trace of A* to vanish, one obtains that 
i Pnmd) = — (^4* , I — Pprod) ■ Hence, one of the separa- 



(A* , Pprod) 

ble states a' = P, 



prod 7 



jyzi(/ - Pprod) violates the 



condition (|13j), which gives the expected contradiction. 
The above theorem leads immediately to the following 



Theorem 2 The measure n(S sep ) of separable states is 
a nonzero one. In particular there exists always some 
e > such that the following inequality holds 



e"- 1 > 0. 



(16) 



As an illustration of the above theorem, let us consider 
the 2 x 2 or 2 x 3 cases (N = 4, 6) for which separability is 
equivalent to the positivity of the partial transposition. 
It is easy to see that the spectrum of the partially trans- 
posed density matrix must belong to the interval [— i , 1] . 



Hence any state of the form g = (1 — p)-L + pg, for an 
arbitrary g and p < 2/(2 + N), has a positively defined 
partial transposition, and thus is separable for the con- 
sidered cases. As the maximal value of p is 1/3, or 1/4, 
it means that the value of e in the above theorem can be 
estimated just by ^, or i for N — 4, 6, respectively. In 
the next section we shall show that those bounds can be 
significantly improved. 



C. Purity and separability 

As we have shown, all states in the small enough neigh- 
borhood of the totally mixed state gi = I/N are sepa- 
rable. On the other hand, we know that in the sub- 
space of all pure states, the measure of separable states 
is equal to zero |2j. It is, therefore, interesting to inves- 
tigate the relationship between entanglement and mix- 
ture of quantum states. A qualitative characterization of 
the degree of mixture is provided by the von Neumann 
entropy Hi(g) = — Tr(glng). Another quantity, called 
participation ratio: 



R(g) 



W) 



(17) 



is often more convenient for calculations. It varies from 
the unity (for pure states) to N (the totally mixed state 
qi) and may be interpreted as an effective number of 
states in the mixture. This quantity, applied in solid 
state physics long time ago |0|, is related to the von 
Neumann-Renyi entropy of order two, H2(g) — lnR(g). 
The latter, called also purity of the state, together with 
other quantum Renyi entropies H q (g) = (ln[Trg 9 ])/(l— q) 
is used, for q =/= 1, as a measure of how much a given 
state is mixed. It has also been applied for the deriva- 
tion of some necessary conditions of separability in Ref. 
]l6[ . In Section VI we shall demonstrate, using numerical 
simulations, that the participation ratio (and other von 
Neumann-Renyi entropies) allows to establish a dualism 
between purity and separability of the states of compos- 
ite systems. In this subsection we use it to calculate a 
natural lower bound on the volume of separable states 
for dimensions N = 4, 6. 

For this purpose consider the following lemma: 



Lemma 3 If the state g satisfies 
R(q) >N-1 



(18) 



where N is the dimension of Tt, then g T2 is positive de- 
fined, i.e. its spectrum s(g T2 ) belongs to the simplex A. 

Proof. - Let us denote by Bn(t, P) the ball in the space 
R N with radius r and center P and by Spf(r, P) its sur- 
face. The condition ( |l8|) is invariant with respect to the 
partial transposition, because Ti(g 2 ) = Tr((g T2 ) 2 ). That 
implies that s(g T2 ) £ Bn(t,z\) with r = -^==,zi = 

(■ir, -h). Let define the (N — l)-dimensional linear 
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manifold M n- i = {x = (xi, ... , ajjv), 2j=i x i — !}■ We 
only need to show that its intersection with the ball is 
included in the simplex A, i.e. that the new (N — 1)- 
dimensional ball B' N1 = Bn{t,z\) n Mn-i C A. This 
can be seen in the following way. It follows from the 
high symmetry of the sphere and the invariance of the 
simplex under cyclic permutations of coordinates, that 
the center of this intersection is again zi. Hence the 
radius r' of Bn—i can be calculated immediately by tak- 
ing the distance from an arbitrary point from the sur- 
face Sjsr(r, Zi) n JAn-i (say, for example from the point 
(0) lv3T> ■■■) jv-i ) to the point zi). It is elementary to 
show that r' — 1/ */ N(N — 1). On the other hand, one 
can calculate the maximal radius r" of the ball of the 
type B'p f _ 1 (r" , Zi) included in A by calculating the min- 
imal distance of zj to the boundary of A. To this aim 
we have to minimize (r") 2 = Yli=i( x i — 1/N) 2 with the 
constraints Yli=i Xj — 0, and x = 0. Using Lagrange 
multipliers we obtain immediately r" — r' , and hence 
Bm-i belongs to A, which ends the proof. 

Now using the explicit expressions for the volume of a 
(N — l)-dimensional ball (V N (r) = 7r j T Li r Jv ^ 1 /r(^-)), 
and for the volume of the simplex A belonging to the 
manifold Mn-i (Va = VN/{N — 1)!), one can obtain 
the lower bound of the volume of states with positive 
partial transposition, 



TN 



(N - 1)!tT 



N^(N-l) — 



T(^±i) 



(19) 



Recalling that for Hilbert spaces of dimensions N = 2x2 
and TV = 2 x 3 the states with positive partial transpo- 
sition are the separable states, Eq. ( |l9| ) leads directly to 
the following theorem: 

Theorem 3 // the participation ratio satisfies R(g) > 
3 (R(g) > 5) for N = 4 (N = 6) then the state g is 
separable. 

Therefore, the measure fi(S sep ) of separable states is re- 
stricted from below by the inequality fj,(S se p) > ^= ~ 



0.302 for TV = 4 and 



625V5 



0.056 for N = 6. 



IV. UPPER BOUNDS ON THE VOLUME OF 
SEPARABLE STATES 

In this section we seek for upper bounds on the volume 
of V sep , or, equivalently, lower bounds on the set of in- 
separable states Vi nse p. Several necessary conditions for 
separability have been recently established with the aid 
of positive maps. We should use them to determine an 
upper bound on V sep . As we shall see, these conditions 
are in some way complementary and can be combined to 
obtain a better estimate of the upper bound of V sep . 
Our first estimate relies on the positivity of the par- 
tial transposition. It is valid for any dimension, but 



we shall apply it to composite systems of dimension 
2x2. Note that if a state has a partial transposition 
which is not positively defined, then the state is neces- 
sarily inseparable. Before proceeding further we should 
first recall the Schmidt decomposition of a pure state 
|*) € H = Hi <g> H 2 , dimWi = N u dimW 2 = N 2 , 
Ni x N 2 = N, 



I*) 



min(Afi,Af 2 ) 

E 

i=l 



(20) 



where |e») ® \fi) form a bi-orthogonal basis (e.j|ej) = 
(fi\fj) — Sij, and < ai < 1 denote the coefficients of the 
Schmidt decomposition with the condition J^. af = 1. 
It is straightforward to see that Pi 2 = (|*)(*|) T2 has 
eigenvalues af for i = 1, . . . , min(Ni,N2) and ztaiaj for 
i ^ 3- 

We can now state our first lemma: 

Lemma 4 If in the range of a state g there exists 
such that 



a = (^ig- 1 !*)- 1 > 



then g is inseparable. 



1 + max^j (a,aj) ' 



(21) 



Proof.- According to Ref. [jl7j, any state g can be ex- 
pressed as: 



g = APq, + (1 - A)g, 



(22) 



where P^ is a projector onto |^) and g is a (positively 
defined) state. Thus: 



g T - - AP^ + (1 - A)g 



'A 



(23) 



Recall that for any iVi, N 2 , the eigenvalues of g T2 belong 
to the interval [—1/2,1]. Let \^ neg ) denote the eigen- 
vector corresponding to the minimal eigenvalue of P^ 2 : 
—maxi^j^aiaj). We have thus 

(^wegl^l^nea) < -A(maxi^(o,-Oj)) + 1 - A < 0, (24) 

because {^ neg \g T2 \^ neg ) < 1. The above inequality im- 
plies that g T2 is not positively defined when condition 
( pl| ) holds, and therefore g is not separable. Note, that 
the lemma (^) can be in particular applied to the eigen- 
vectors of g: 

Lemma 5 If g has an eigenvector corresponding to 
the eigenvalue A such that the condition (jITp holds, then 
g is inseparable. 

The eigenvalue A can fulfill the above condition if and 
only if it is the largest eigenvalue, because it must be 
larger than 2/3. The corresponding normalized eigen- 
vector, however, is absolutely arbitrary, and according 
to invariant measure on the group it can be generated 
simply by a uniform probability on the N dimensional 
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unit sphere. This implies, as we shall see below, that 
the Schmidt coefficients a* are also absolutely arbitrary 
and distributed uniformly on the octant of the \/~N- 
dimensional sphere. 

Consider the Ni = N 2 = K case. Any vector in 
N = K 2 -dimensional space from the unit sphere can 
be represented by a row of complex numbers Xi with the 
condition J^. = 1. In any product basis, we can 
view it as a K x K matrix Cy with i,j — 1, . . . , K, and 
with the condition Tr(C^C)) = 1. We seek the uniform 
distribution on the set of such matrices. But from the 
polar decomposition theorem any matrix of such a type 
can be represented in the form: 



C = U'DU 



(25) 



where U' and U are some unitary matrices, while D is di- 
agonal matrix with nonnegative elements (eigenvalues). 
These eigenvalues are nothing else but at. The reason is 
that the above form, which is the analog of the spectral 
decomposition of the hermitian matrix, is at the same 
time the Schmidt decomposition written in the matrix 
notation. In our case (taking into account the above 
mentioned trace condition), the spectrum of D is repre- 
sented by the point belonging to the octant area of the 
sphere. These leads to the measure 



//(*) = v(U\K)) x u(U(K)) x n(D), 



(26) 



where the first two measures are Haar measures on the 
unitary group U(K), and the last one is the uniform 
(Lebesgue) measure on the octant of the ball in K di- 
mensional space. Similar results can be straightforwardly 
generalized for the cases N\ ^ N 2 . 

If one calculates now the measure ( prf ) for |^), and 
combines it with the uniform measure on the simplex , 
one could estimate an upper bound of separable states: 

H{S S ep) < 1 ~ 

J e(maxAi - (1 + max i ^ J (a i a i )) _1 )d^'(*)d^A, (27) 

where denotes the Heavyside function. The double 
integration over the unitary groups that contains 
can be easily performed, since neither Aj nor a, depend 
on the direction of |\&). This is the first qualitative ar- 
gument that the measure of inseparable states does not 
vanish. 

Moreover, recently a new separability condition 
have been introduced with the aid of positive maps con- 
dition |To|]: if the state g is separable then I ® gi — g 
must be positive, where g\ is the reduced density matrix. 
It implies for any |*) that Tr((I ® Qi)Pn) > Tr(gP^). 
Straightforward estimation tells us that for any separa- 
ble state it must hold that < max of, where di 

i 

are again the Schmidt decomposition coefficients of ^ . 
That implies a lemma analogous to lemma (0). 



Lemma 6 If in the range of a state g there exists \9) 
such that 



A = (*|e|*> > max of 



(28) 



then g is inseparable. 



This lemma is neither stronger nor weaker than the 
lemma If we apply it to eigenvectors of g, however, 
the relevant eigenvalue need not be the maximal. In the 
case 2 x 2 we can combine both conditions (Lemma ||-|(]) 
to obtain a better estimate on the upper bound of /i(Sep) 



1 - [i{S sep ) > 



V A V C 



oct Jo 



dAi 



l-Ai 



dkn 



I-A1-A2 



/ dai I da 2 9(Ai - (1 

Jai>0 Ja 2 >0 



dA 3 



0(Ai — max(aj, a^))^ ! + a 



1) 



(29) 



Notice that in the above expression the first three inte- 
grals are over the eigenvalues of g that are located in the 
simplex A, whereas the remaining two integrals are on 
the eigenvalues of D from the octant area of the sphere. 
The integrals can be calculated analytically but the re- 
sulting expressions are very complex. After a tedious, 
but straightforward calculation we obtain: 



KSsep) < 0.863 



(30) 



In general for an arbitrary dimension N = Nix N 2 and 
K = min(A r i, A/^); 1 — [i(S sep ) can be estimated from 
below by a bound b(Ni,N 2 ) using the above method. 
This bound, on the other hand, can be estimated from 
above by the volume of the "corners" of the simplex A of 
the sides 1 — regardless the uniform measure on pure 
states by putting it equal to unity there. This follows 
from the fact that the condition (P8|) can only then be 
fulfilled, when an eigenvalue of g is larger than max a 2 > 

i 

1/K. The relative volume of such corners equals N(l — 
j^) 1 ^^ 1 - Keeping in mind the formula ( pfj| ) the above 
simple estimation leads to the corollary: 

Corollary 1 Consider a quantum system G Hi ® Ti 2 , 
where dimHi = Ni, dimH 2 = N%, N = Ni x N 2 and 
K = min(Afi, AT 2 ). Then using Lemae 4~6 the volume of 
separable states is restricted from above by: 



KSsep) <l~b(Ni,N 2 ) 



where: 



b(Ni,N 2 )<(l--) 



N-l 



(31) 
(32) 



The volume (|3^), however, converges asymptotically 
to the value 1 as Ni,N 2 grow, so that in the limit of 
large N we get a trivial result /x(5 sep ) < 1. At the same 
time, the numerical results which we shall present subse- 
quently strongly suggest that there should exist an upper 
bound for ^,(S sep ) converging to zero. So far, the rigorous 
proof that in the limit of higher dimensions fi(S sep ) — ► 
remains an open problem. 
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V. INSEPARABLE STATES WITH POSITIVE 
PARTIAL TRANSPOSITION 



VI. NUMERICAL RESULTS 



As it was mentioned in the introduction for N > 8 
there are states which are inseparable but have posi- 
tive partial transposition ||,[l(|. Moreover, it has been 
recently shown that all states of such type represent 
"bound" entanglement in the sense that they cannot be 
distilled to the singlet form M . The immediate question 
that arises is how frequently such peculiar states appear 
in the set of all the states of a given composite system. 
This question is related to the role of time reversal in the 
context of entanglement of mixed states |],|7j Below we 
provide a qualitative argument that the volume of the set 
of those states is also nonzero: 

Lemma 7 For N > 8 the set of inseparable states with 
positive partial transposition includes a nonempty ball. 

Proof.- Consider the two sets of quantum states for 
some composite system: the set of separable states S se p 
and the set of states with positive partial transposition 
T. The first of them is convex and compact. The second 
one is convex set [| Since positivity of partial trans- 
position is necessary for separability, we have obviously 
Ssep € T. Consider any state a belonging to T but not 
to S sep (we know that for N > 8 such states exist). Let 
us take the ball B(r,Qj) around the maximally chaotic 
state Qi such that the whole B(r, gj) belongs to S se p (the 
ball can be in principle defined in an arbitrary norm as 
all norms are topologically equivalent, c.f. section III). 
Obviously such a ball exists; otherwise gi would belong 
to the boundary of S se p ■ It can be shown that the latter 
would contradict the fact of a nonzero volume (see sec- 
tion III). Consider the sequence of balls obtained from 
B{r, gi) by a translation of the center r, and a rescaling: 
B n = B(Z-,(l—^)a+~Qi). Some of B„'s have to include 
no separable states. Otherwise, if any B n included some 
separable states e.g. gff p , by virtue of compactivity of 
Ssep, the state a would be separable as the limit of a se- 
quence of separable states g s ^ v . Hence, some B no do not 
belong to set of separable states. But on the other hand, 
any state from B no is a convex combination of elements 
of T. Thus the whole ball B no belongs to T, as the latter 
is a convex set. In this way we have shown that there is 
some ball B„ C T and does not intersect with S se p- But 
the inseparable states with positive partial transposition 
are just the ones belonging to T and not to S sep - This 
ends the proof of the lemma. 



3 In fact it is also compact but we shall not need this property 
here. 



In this section we provide rather precise estimates of 
the measure n(S se p) for N = 4, 6, as well as upper bounds 
of such measure for N > 6. Our results are obtained nu- 
merically. This section is self-contained, in the sense that 
can be read directly without passing though the more 
technical sections III-V. It also includes the results ob- 
tained in these previous sections. 

A. Estimation of volume of separable states 

Our goal is to estimate the volume of the set of sep- 
arable states fi(Ssep)- For simplicity we discuss states 
consisting of 2 subsystems. Any state describing a mix- 
ture of N\ and ^-dimensional subspaces may be repre- 
sented by a positive defined (NiN 2 ) x (N1N2) hermitian 
matrix g with trace equal to unity. g T2 denotes, as be- 
fore, the matrix partially transposed with respect to the 
second subsystem. As mentioned previously, if g is sepa- 
rable, then necessarily g T2 is positive plf . Moreover, for 
the simplest 2x2 and 2x3 problems this condition is 
also a sufficient one jl(J , which is not true in the general 
case N > 8 Therefore, the set of separable states 

is a subset of states with positive partial transpositions. 
Thus, in order to estimate the volume of separable states 
from above, it is sufficient to find the volume of the set 
of states with positive partial transpositions. 

Let us remind that according to the section II any state 
(any density matrix) can be represented in a family V of 
complete sets of orthogonal projectors, and the simplex 
A representing all possible spectra, 

S = V x A. (33) 

On the other hand, any element of V can be represented 
by a unitary transformation, and any element of A, as a 
diagonal matrix D with the matrix elements Ay = (JyAj, 

such that they fulfill ^—^Ai = I. Such representation 
corresponds to the form: 

g=UDU f . (34) 

Thus, a uniform distribution on the set of all density 
matrices represented by (^3|) is constructed naturally by 
postulating an uniform distribution on unitary transfor- 
mations U(N) (the Haar measure), and an uniform dis- 
tribution on diagonal matrices D. 

We have calculated numerically the volume of the set 
of matrices with positive partial transpose, and we have 
estimated in such a way the volume of the separable 
states. An algorithm to generate random U(N) matrices 
was recently given in Refs. The random diag- 

onal matrix D fulfills that Yli=i Aj = 1, so the vector 
A = (Ai, . . . , Ajy) is localized on the N — 1 dimensional 
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simplex A (see section II). Physically speaking, no com- 
ponent of this vector is distinguished in any sense. Ran- 
dom vectors A are thus generated uniformly on this sub- 
space according to the simple method described in the 
Appendix A. 

The numeric algorithm is then straightforward: firstly 
we generate random density matrices of any size N — 
Ni x N 2l secondly we construct their partial transpo- 
sitions and finally we diagonalize them and we check 
whether their eigenvalues Xf, i = 1, . . . , N are all positive. 
This procedure has been repeated several millions times 
in order to obtain the accuracy of the order of 1/1000. 

We should recall that in previous sections we have 
obtained rigorous analytical lower and upper bounds of 
n(Ssep), i-e. we have proven that < fi(S sep ) < 1. The 
lower bounds follow from the fact that the states suf- 
ficiently close to the totally mixed state are separable 
for all N. We have also shown that states which have 
sufficiently large participation ratio (i.e. which are suffi- 
ciently mixed) have a positive partial transpose, and are 
thus separable for N = 4, 6. On the other hand, the up- 
per bounds come from the fact that matrices with large 
eigenvalues corresponding to an entangled eigenvector, 
are necessarily inseparable. Let us denote the measure 
li{S sep ) in the iVi x N 2 case by -P/v lX Af 2 - Our analytical 
bounds for the cases 2x2 and 2x3 are summarize below: 

0.302 < P 2x2 < 0.863, 

0.056 < P 2x3 . (35) 

Our numerical results agree with these bounds but to our 
surprise the probability that a mixed state g <E H 2 x H 2 
is separable exceeds 50%. The results are 

P 2x2 « 0.632 ±0.002 and P 2x3 w 0.384 ± 0.002. (36) 

For higher dimensions our results are summarized in Fig- 
ure 1. This figure displays the probability Pjy that the 
partially transposed matrix g T2 is positive as a function 
of N = Ni x N 2 . For N = 4 and N = 6 this is just the 
required probability of encountering a separable state, 
while for N > 6 it gives an upper bound for this quan- 
tity. 

Due to symmetry of the problem Pj^ 1 x jy 2 must be equal 
to Pjv 2 xJVi- Numerical results strongly suggest that this 
quantity depends only on the product N — Ni x N 2 , 
eg. -P 2 x6 — p3x4- Moreover, the dependence can be 
well reproduced by an exponential decay. The best fit 
gives Pm ~ 1.8e~ a26JV . We conjecture, therefore, that 
the measure of separable states decreases exponentially 
with the size of the system in consideration. 

B. Purity versus separability 

Here we would like to illustrate the physical connection 
between the participation ratio and entanglement, which 
was already discussed in Section III.C. Recalling that the 



participation ratio: R(g) = 1/Tr(g 2 ), gives a characteri- 
zation of the degree of mixture and can be interpreted as 
the effective number of states on the mixture. We have 
demonstrated that if the state g has a sufficiently large 
participation ratio, or equivalently a sufficiently low von 
Neuman-Renyi entropy H 2 (g) = IvlR(q), then it partial 
tranposed is always positive. This holds for any arbitrary 
N, so in particular it means it is separable for N = 4, 6. 
A more precise estimate can be performed numerically. 

For example consider the N = 4-dimensional Hilbcrt 
space. A manifold of the constant participation ratio R 
is given by the ellipsoid in the space of eigenvalues: 

A? + A 2 + A 2 + (1 - Ai - A 2 - A3) 2 = 1/R. 

The probability distribution P(R) obtained numerically 
using the natural uniform measure on the 3-dimensional 
simplex is plotted in Fig 2(a). It corresponds to the rel- 
ative volume of the cross-section of a 3D hypersphere of 
radius Rr 1 / 2 centered at (0,0,0,0) with the simplex de- 
fined by a condition, Ai + A 2 + A 3 + A 4 = 1. For N = 4 
and R > 3 we obtain P(R) = 6nR- 2 ^/l/R - 1/4. 

We have generated randomly a million points in the 
3D simplex, computed the corresponding participation 
ratio, rotated the corresponding state by a random uni- 
tary matrix U and checked, whether the generated state 
is separable. This procedure allows us to investigate the 
dependence of the probability of the separable states on 
the participation ratio. Our numerical results are sum- 
marized in Fig 2 (b), and again, are fully compatible 
with the theorems and lemmas obtained in the previous 
sections. 

Similar results and compatibility are obtained for the 
case N = 6. In this case we deal with a 5D simplex. Nu- 
merical data, displayed in Fig. 3, support the general fact 
that the quantum states with R(g) > N — 1 have posi- 
tive partial transpose (i.e. are separable for N — 4,6) in 
the A^-dimensional Hilbert space. It is interesting to note 
that the relative amount of separable states unambigu- 
ously increases as the participation ratio grows. More- 
over, the mean 'degree of entanglement' (t), defined in 
the Appendix B, decreases monotonically with R. This 
illustrates that for composite systems there exist a dual- 
ism between the two quantities: purity and separability 
of the state. The purest a quantum state is, the small- 
est its probability of being separable. This conclusion is 
supported by Fig. 4 where we plot the dependence of the 
probability of finding a separable state P sep on the values 
of several quantum Renyi entropies H q (g). For 3 / 1 the 
Renyi entropy is defined as: H q (g) := (ln[Trg 9 ])/(l — q), 
while in the limit q — > 1 it tends to the standard von Neu- 
mann entropy H\{q) = — Trg In g. We have not been able 
to generalize the rigorous the results for q = 2, for which 
large entropy implies necessarily positivity of the partial 
transpose. Nevertheless, the numerical results suggest 
that similar result holds for arbitrary values of the Renyi 
parameter q. All states with sufficiently large H q {g) are 
separable as shown in figure 4a. Can this fact be used to 
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obtain a better lower bound for -P/v iX jv 2 ? Fig. 4b sug- 
gests that it is not the case: the cumulative probability 
P of H q (g) attains more or less the same value for all q 
at the point at which P sep becomes ~ 1. 

VII. CONCLUSIONS 

Summarizing, we have developed a measure theoretical 
approach to the separability-inseparability problem. To 
this aim, we have proposed a natural measure in the space 
of density matrices g on the iV-dimensional space. We 
have proven that, under this measure, the set of separable 
states has a nonzero volume although this volume is not 
maximal in the set of all states. Analytical lower and 
upper bounds of this volume have been found for N — 
2x2, and N — 2 x 3 cases. We have also provided 
qualitative evidence that for N > 8 the peculiar set of 
inseparable states with positive partial transposition has, 
under this measure, a nonzero volume. 

We have used Monte Carlo simulations to estimate 
with much higher precision the volume of separable 
states. Our numerical simulations give strong evidence 
that this volume decreases exponentially with the dimen- 
sion of the composite quantum system. Finally, we have 
also discussed the dualism between purity and separabil- 
ity, and have shown that while entanglement is typical of 
pure states, separability is rather connected with quan- 
tum mixtures. 

Several questions concerning this subject remain still 
as open problems and so far, we have not been able to 
prove them rigorously. Particularly challenging are the 
two following related questions: 

• Does the volume of the set of separable states goes 
really to zero as the dimension of the composite 
system N grows, and how fast? 

• Has the set of separable states really a volume 
strictly smaller than the volume of the set of states 
with a positive partial transpose? 
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APPENDIX A: GENERATION OF UNIFORM 
DISTRIBUTION ON THE SIMPLEX 

Our aim is to construct the uniform distribution of the 
points on the manifold given by At = 1, where each 
component A, is non negative. Let = 1, . . . , N — 1 
be independent random numbers generated uniformly in 
the interval (0, 1). We start with a uniform distributions 
inside the N — 1 dimensional simplex A„_i defined by 
X3 i= ^ Aj < 1. Its volume is proportional to the product 

Ylk=i x^~ 1 ~ k dxk = Ylk=i d[x^~ k ], what enables us to 
find the required densities for each component. Since 
the vertex of the simplex Ajv-i is situated at {0, . . . , 0}, 
the largest weight corresponds to the small values of X\ . 
Therefore: 

Ai = i-er\ 

A 2 = [i-e 2 wts ](i-A 1 ), 
A fc = [i-er*](i-X>)> 

i=l 



AT-2 
i=l 

Eventually the last component Ajv is already determined 
as: 

JV-l 

A N = I - A,. 

i=l 

The vector A = {Ai,...,Aj\r} constructed in this way 
is distributed uniformly in the requested subspace. An 
alternative procedure, albeit more time consuming, is to 
take any vector of an N x N auxiliary random unitary 
matrix V and obtain the random vector as A& = |Vfej| 2 
with arbitrary j. 

APPENDIX B: THE AVERAGED "DEGREE OF 
ENTANGLEMENT" 

The problem of defining a quantity capable to measure 
a "degree of entanglement" is a subject of several recent 
studies pj|-|22]|. Let us define for a given density matrix 
g the quantity 

N 

i=i 
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where A^, i = 1, . . . , N denote the eigenvalues of the par- 
tially transposed matrix g T2 . For any separable matrix 
all eigenvalues are positive, its trace is equal to unity and 
t equals to zero. On the other hand, for the maximally 
entangled states belonging to 2 x 2 system the spectrum 
of eigenvalues A' consists of {—1/2, 1/2, 1/2, 1/2}, so that 
t = 1. Moreover, for the often studied 2x2 Werner 
states ||] depending on the parameter x, the quantity 
t vanishes for x < 2/3 (separable states) and equals to 
t = (3x - 2)/(4-3x) for entangled states (2/3 < x < 1). 

We could not resist the temptation to investigate the 
mean value of t averaged over random density matrices 
generated as described above. For the 2x2 problem the 
mean value (t) equals to 0.057 and increases to 0.076 for 
the 2x3 problem. For large systems this quantity seems 
to saturate at t ~ 0.10, as the ratio of the matrices with 
positive values of t (some eigenvalues of g T2 are negative) 
tends to unity. Moreover, as shown in Fig. 2b and Fig 3b, 
the average degree of entanglement (t) decreases mono- 
tonically with the participation ratio R, what provides 
a quantitative characterization of the relation between 
entanglement and purity of mixed quantum states. 
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FIG. 2. Purity and separability in iV=4-dimensional 
Hilbert space: a) probability of finding quantum state with 
a participation ratio R; b) probability of finding a separable 
state P2x2 as a function of the participation ratio R (crosses). 
All states beyond the dashed vertical line placed at R — N — 1 
are separable. Circles show the mean entanglement (t) , as de- 
fined in Appendix B. 



FIG. 3. Same as in Fig. 2 for N = 6. 



FIG. 4. (a) Probability of separable states P2x2 for N = 4 
vs von Neumann- Renyi entropies H q 6 [0, ln(4)] for q = 1(*), 
2(o), 3(A), and 10(o). (b) Integrated distribution function 
D(H q ). Vertical line drawn at D ~ 0.7 corresponds to these 
values of H q for which P sep = 1. 
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